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ABSTRACT 

Elemental abundance effects in active coronae have eluded our understanding for almost three decades, since the discovery of the first 
ionization potential (FIP) effect on the sun. The goal of this paper is to monitor the same coronal structures over a time interval of 
six days and resolve active regions on a stellar corona through rotational modulation. We report on four iso-phase X-ray spectroscopic 
observations of the RS CVn binary El Eri with XMM-Newton, carried out approximately every two days, to match the rotation period 
of El Eri. We present an analysis of the thermal and chemical structure of the El Eri corona as it evolves over the six days. Although 
the corona is rather steady in its temperature distribution, the emission measure and FIP bias both vary and seem to be correlated. An 
active region, predating the beginning of the campaign, repeatedly enters into our view at the same phase as it rotates from beyond 
the stellar limb. As a result, the abundances tend slightly, but consistently, to increase for high FIP elements (an inverse FIP effect) 
with phase. We estimate the abundance increase of high FIP elements in the active region to be of about 75% over the coronal mean. 
This observed fractionation of elements in an active region on time scales of days provides circumstantial clues regarding the element 
enrichment mechanism of non-flaring stellar coronae. 
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1. Introduction 

The best studied cool-star corona is that of our own sun, which 
also provides the only spatially resolved stellar coronal source 
of extreme ultraviolet (EUV) and X-rays. From X-ray images 
of the sun we know that the solar corona in its active state is 
non-uniform and that a significant fraction of the emission in 
these bands comes from magnetic loop structures concentrated 
in active regions around photospheric sun-spots. Evidently, these 
structures evolve in brightness on time scales of hours. The mor- 
phology as well as the metal abundances vary on longer time 
scales of days. Coronal structures on the Sun may survive for 
a few weeks, although tracking them continuously from the 
ground is difficult due to the ~24.5 day solar rotation period. 

Active coronae of other cool stars can be up to five orders 
of magnitude more luminous and an order of magnitude hot- 
ter than the solar corona. Although unresolved, they show many 
characteristics similar to those of the solar corona, and are there- 
fore believed to be structured and to evolve in a similar way. For 
the most part, late-type stars have envelopes in which convec- 
tion and differential rotation interact to drive a dynamo, which 
produces magnetic flux tubes. The flux tubes float into the up- 
per atmosphere and eventually the corona. It is generally un- 
clear, however, whether stellar coronae actually have larger cov- 
ering of their stellar photosphere, or even qualitatively different 
morphologies, such as inter-binary magnetic structures. Stellar 
coronae show high variability, often attributed to continuous 


micro-flaring, which is observed on the sun, but perhaps with 
different statistics (e.g., Caramazza et al. 2007). 

Solar coronal chemical abundances are different from the 
solar photospheric composition. The solar corona itself is non- 
uniform in that respect. While in coronal holes, the abundances 
are close to photospheric, abundances in active regions and in 
the solar wind are significantly different. The abundance pat- 
terns depend on the first ionization potential (FIP; Meyer 1985; 
Feldman 1992; Sheeley 1995). Elements with FIP < 10 eV 
(low FIP) are enriched in the corona relative to elements with 
FIP > 12 eV (high FIP). Averaging over the entire solar disk, 
the enrichment factor is about 4 (Feldman & Widing 2002, and 
references therein). Such a FIP bias, in varying intensities, has 
been reported for many coronal structures, from closed magnetic 
loops to streamers (see reviews by Kohl et al. 2006; Doschek & 
Feldman 2010). Currently there is no generally accepted model 
that explains the FIP bias. 

The FIP effect became even more puzzling when high reso- 
lution X-ray spectra from XMM-Newton and Chandra revealed 
that active stellar coronae do not follow the solar FIP pattern. In 
some cases, the high-FIP elements are even enriched compared 
to the low-FIP ones, an effect labeled the inverse-FIP (IFIP) ef- 
fect (Brinkman et al. 2001). Eater studies revealed that FIP and 
IFIP biases are correlated with coronal activity (and age): Highly 
active stars show an IFIP effect, while less active coronae have a 
solar FIP bias (Audard et al. 2003; Telleschi et al. 2005), imply- 
ing a transition on stellar-evolution time scales of Gyr. 
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Table 1 . Observing log. 


Obs. ID 

Start time 

End time 

Exposure [ks] 

060317201 

2009-08-0 1T00:02:50 

2009-08-0 1T09: 13:05 

32.84 

060317301 

2009-08-02T23:06: 16 

2009-08-03T09:24:46 

36.92 

060317401 

2009-08-04T22:20:49 

2009-08-05T07 :42 :4 1 

33.54 

060317501 

2009-08-06T20:59:39 

2009-08-07T07 :3 1 : 30 

37.48 


On the face of it, the correlation with activity supported pre- 
vious suggestions that abundances may be linked to flares (e.g., 
Giidel et al. 1999, 2004; Osten et al. 2003; Nordon et al. 2006). 
However, direct systematic measurements of abundance varia- 
tions during flares show that regardless of coronal FIP/IFIP bias, 
during large flares the abundances tend to photospheric values 
(Nordon & Behar 2007, 2008). The mechanism for FIP/IFIP is, 
therefore, related either to the long-term evolution of coronal 
structures, or to the continuous micro-flaring that produces the 
quiescence emission, but not to large flares. 

In recent years, a model has been proposed in which abun- 
dance fractionation is driven by magneto hydrodynamics (MHD) 
waves (Laming 2004, 2009, 2012). This “pondermotive” force is 
controlled by the energy density of MHD waves at the bases of 
the coronal loops. Wood et al. (2012) studied a sample of X-ray 
faint (L x < 10 29 erg s -1 ) dwarfs and found no correlation with 
X-ray luminosity or surface flux, but a good correlation with 
spectral type instead. According to this relation M dwarfs are 
IFIP biased while early G dwarfs are FIP biased. Wood et al. 
(2012) claim that this shift from FIP to IFIP can be explained 
by the inversion of the pondermotive force at certain conditions. 
The latter study is limited to X-ray faint dwarfs and may not be 
applicable to the typical highly active binaries from other stud- 
ies that are more than two order of magnitude more luminous in 
X-rays, emitted mostly by a giant or sub-giant star (RSCVn type 
systems). 

Despite the assumed similarity between solar and stellar 
coronae, we have lacked a study of the evolution of a stellar 
corona on time scales of days, which is the relevant time for 
coronal loop evolution on the sun. The goal of the observations 
presented in this paper is to seek and measure the variability of 
temperatures and abundances in an active corona over several 
days. In order to avoid confusion between rotational variability, 
which is inevitable due to the rapid rotation typical of active stel- 
lar targets and the corona being composed of discrete active re- 
gions as on the sun, and genuine coronal-activity variability, the 
present observing campaign is designed to return to the target 
each time at exactly the same orbital phase. 

1.1. El Eh 

To carry out a few iso-phase observations, we chose the tar- 
get of El Eridani (El Eri, HD 26337), an X-ray bright (L x = 
10 31 ' 1 erg/s) RS CVn binary, whose orbital period is 1.95 days 
(168.5 ks), which is very similar to the exactly 2.0 day orbital pe- 
riod of XMM-Newton around Earth. This enabled XMM-Newton 
to observe the exact same orbital phase of El Eri during four 
consecutive orbits, and for sufficiently long exposures each time 
( — 35 ks). El Eri is a nearby (50.8 pc, van Leeuwen 2007), rapidly 
rotating (psin i = 51 kms -1 ), non-eclipsing, rotationally locked, 
close RS CVn binary with an inclination angle to the rotation 
axis of ~56° ± 4.5° (Washuettl et al. 2009a). These properties 
make El Eri a popular target for Doppler imaging and starspot 
monitoring that indicate the primary star to be highly active. 
El Eri consists of a solar-mass G5IV primary and a much fainter 


M4-5 secondary, with a possible 19-year orbit of a third body 
(Washuettl et al. 2009a). 

Long term photometric monitoring of El Eri seemed to indi- 
cate a ~1 1 yr brightness cycle that was conjectured to be an ana- 
log of the solar activity cycle (Olah & Strassmeier 2002), but un- 
like the sun, long term Doppler imaging of El Eri by Washuettl 
et al. (2009b) shows little variability of its cool starspots 
over 1 1 yr. As for the hot corona, Osten et al. (2002) observed 
El Eri with ASCA (and EIJVE) and fitted its spectrum with a 
2T model of kT = 0.72 keV and 2.2 keV, but the abundances 
could not be tightly constrained with ASCA. The XMM-Newton 
CCD spectra of a short flare on El Eri on 2009-08-06 during the 
present campaign was analyzed by Pandey & Singh (2012), who 
suggested the flare originated from a loop the size of ~R, (pri- 
mary), which is about a tenth of the binary separation distance. 

2. Observations, data, and analysis 

In Table 1 we list the log of the four XMM-Newton observations 
of El Eri. Note the start times for each observation, which are 
separated by approximately the El Eri orbital period of 168.5 ks. 

2.1. Data reduction 

XMM-Newton data were reduced using the Science Analysis 
Software (SAS - ver. 10.0.0) in the standard processing chains 
as described in the ABC Guide to XMM-Newton Data Analysis 1 , 
and using the most updated calibration files (date 22.06.2010). 
Pileup was checked using the epatplot routine of the SAS soft- 
ware package. No significant pileup was detected in either of the 
EPIC instruments and thus events from single to quadruple were 
retained. Background was low during all the observations and no 
significant amount of time was lost. 

2.2. Spectral analysis 

We measure the temperature and abundance properties of the 
El Eri corona by fitting thermal plasma models to the observed 
spectra. XMM-Newton comprises several different instruments. 
The longer wavelengths (6-38 A) are covered by the two high 
resolution reflection grating spectrometers (RGS) that resolve 
the emission lines. The high energy end up to 15 keV is covered 
by the two EPIC-MOS and single EPIC-pn CCDs, that provide 
lower spectral resolution. 

To study the thermal and chemical composition of the source 
plasma, it is imperative to first properly account for its emission 
measure distribution (EMD) that is defined as 

EMD(T) = n e n H dV/dT. (1) 

where n e and « H are the electron and proton number densities, 
respectively, and dV/dT is the temperature distribution of the 
plasma volume. The EMD is usually characterized by several 

1 See http://heasarc . gsfc.nasa.gov/docs/xmm/abc 
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discrete temperature bins. Here, we use seven temperature bins 
with identical abundances. To allow for a uniform analysis of 
the four observations, the temperatures of the middle five com- 
ponents kT = 0.7, 1.0, 1.5, 2.0, and 2.5 keV are held fixed, 
while the lowest and highest temperatures are allowed to vary 
and are determined by fitting the model to the observed spec- 
tra. We fit the observed spectra using the XSPEC (Arnaud 1996) 
software package and a VAPEC (Smith et al. 2001) plasma spec- 
tral model. For high statistics, we fit the two RGS, EPIC-pn, and 
EPIC-MOS1 spectra of each observation (orbit) simultaneously. 
EPIC-MOS2 was used in timing mode and its spectra are sys- 
tematically offset by a small amount from EPIC-pn and EPIC- 
MOS1 that agree with each other very well. For this reason, we 
exclude it from the spectral fitting. We ignore the very low en- 
ergy ( E < 0.35 keV) portions of the CCD spectra. Visual inspec- 
tion of each RGS spectrum indicates the models do an overall 
good job of fitting the lines. See the appendix for the spectra and 
fitted model. The absorption column density towards El Eri is 
negligible (Ah = 1.2 x 10 20 cm -2 , Pandey & Singh 2012) and 
was not included in the model. 

The seven plasma components provide only a general idea of 
the temperature distribution in the plasma. It is known (see most 
recently, Testa et al. 2012) that the EMD reconstruction should 
be treated with caution. More specifically, the EMD is highly 
degenerate between adjacent components. Therefore, we keep 
track of the covariance matrices (between parameters), which are 
fully taken into account when propagating the errors of all de- 
rived quantities. Due to this degeneracy, the more important and 
reliable quantity is the cumulative, or progressively integrated 
emission measure (J EMD AT which we denote by IEM. The 
calculated errors for the IEM are more meaningful than those 
of the EMD, as they are essentially independent of the choice 
of temperature bins, as opposed to errors on the EMD that are 
affected by the arbitrary choice of temperatures and by the prox- 
imity of the other temperature components. Monte Carlo (MC) 
simulations are used for error estimates. The MC simulation ran- 
domly rolls a new set of EMD values according to the correlation 
matrix and individual-component errors obtained from the best- 
fit solution. The central 68% of the distribution for each point in 
the IEM is taken as the equivalent of the l-<x confidence inter- 
val. For a more detailed description of the error propagation see 
Appendix B. 


3. Results 

3.1. Light curves 

By synchronizing the XMM-Newton exposures with the rotation 
phase of El Eri, we are able to observe the same stellar surface 
during each visit, and to thus remove most of the rotation ef- 
fects. The EPIC-pn light curves of the four exposures are plotted 
in Fig. 1 versus the rotaional phase. Only about 20% of the orbit, 
i.e., 35 ks or ~0.4 days, is covered during each visit, and consec- 
utive visits are separated by the remainder of the El Eri period 
(~2 days). 

Between phases 0-0. 1 , all 4 exposures show a steady count 
rate with some variability on time scales of minutes, but at 
a low amplitude that is not much above the Poisson noise. 
This variability is typical of active coronae, and has been as- 
cribed to micro-flaring. The second half of each observation, i.e., 
phases 0. 1-0.2, features a gradual rise in flux that is interpreted 
as the recurrent appearance of an active region on that part of the 
star. 



Phase 

Fig. 1 . EPIC-pn light curves during the four consecutive El Eri rota- 
tion orbits, indicated by sequential ID numbers. The time axis has been 
folded to represent the rotation phase within the 1.95 days period of 
El Eri. 


This brightening with phase is most dramatic during the third 
observation (401), ~4 days into the monitoring, but should still 
not be confused with a coronal flare that would tend to sharply 
rise within a few ks and quickly decay over ~ 10 ks. Indeed, dur- 
ing the last exposure (Obs. 501), a large flare is superimposed 
on the gradually increasing light curve, rising sharply (~1.5 ks) 
about half way through the exposure and decaying exponentially 
on a time scale of ~4.8 ks. For analysis purposes, we therefore 
further split observation 501 into its quiescence and flare peri- 
ods, as indicated in Fig. 1, and refer to them as Obs. 501Q and 
501F, respectively. 

The mean count rate between visits (~ 1 30 ks) also varies, 
increasing within the first 4 days, between Obs. 201 and 401 
by ~20%, and then decreasing by ~30% by the beginning of 
Obs. 501. Interestingly, just before the flare and ~5 days into the 
monitoring, the steady quiescent level is observed to be back to 
its initial level when the monitoring started (Obs. 201, 301). 

3.2. Temperature 

In Fig. 2 we plot the EMD(T) for each observation based on our 
7-temperatures best-fitted models, where each EMD bin is due to 
a single temperature component. The width of each plotted bin 
is arbitrary and designed to give an impression of a continuous 
distribution, although in reality these are iso-T components. As 
with the light curves, the EMD of observations 201, 301, and 
50 IQ are fairly similar, primarily showing emission of plasma 
up to kT — 2 keV, while the EMD of observation 401 features 
some excess high-T emission up to kT - 4 keV. During the 
flare (Obs. 501F), significant emission measure is observed at 
kT > 4 keV. 

The IEMs for each observation are then plotted in Fig. 3. The 
IEM in observations 201, 301 and 50 IQ are consistent with no 
variation in the total emission measure (see the high temperature 
end of Fig. 3), and show very small variations below 1 keV. On 
the other hand, the El Eri corona during observation 401 has less 
emission measure below 1 keV, but a higher total, implying a sig- 
nificant excess of plasma at higher temperatures, mostly between 
2. 5-3. 5 keV. During the short flare (Obs. 501F), one can see an 
excess of emission measure between 1-2 keV and another small 
component above 3 keV that brings the total emission measure to 
a level ~40% higher than the quiescent states of Obs. 201, 301, 
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Fig. 2. Emission measure distribution as obtained by the best-fitted 
7T models of each observation. 



Fig. 3. Cumulatively integrated emission measure distribution IEM cal- 
culated from the EMDs of Fig. 2. Errors take the degeneracy between 
the different temperature components of the models into account. 


and the pre-flare 50 IQ, which is also somewhat lower than dur- 
ing Obs. 401. In short, the emission measure of the quiescent 
corona of El Eri can vary between rotation orbits, and on time 
scale of days, but mostly in the 1-2 keV temperature regime. 
Excess emission measure of hotter plasma up to 4-5 keV is ob- 
served both in active regions (Obs. 401), and more dramatically 
during the flare (Obs. 501F). 

A concise way to quantify the plasma conditions is through 
its mean temperature ( kT ), weighted by the EMD of each com- 
ponent. The errors are again computed using an MC simulation, 
while taking the correlations between the fitted EMDs of the 
various model components into account. The resulting (kT) val- 
ues are plotted in Fig. 4 for each observation and against time 
elapsed from the beginning of the first observation. A similar 
pattern as before emerges. When the source is brighter (e.g., 
Obs. 401) it is also hotter on average. Obs. 201 and 301 do not 
strictly follow this pattern, but the temperature decrease from 
Obs. 301 to Obs. 201 is insignificant. In terms of (kT), we ob- 
serve temperature variations on time scales of days, and most 
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Fig. 4. Variability of the EMD-weighted mean temperature. Horizontal 
error bars indicate the exposure time of each observation. 



FIP [eV] 


Fig. 5. Coronal abundances relative to Fe and expressed relative to the 
solar X/Fe abundance ratio (Asplund et al. 2009). Points from different 
observations are progressively shifted to the right along the horizontal 
axis in the order in which they were taken (~2 days apart) to enhance 
visibility and to emphasize the time variability. 


dramatically before and during the flare. Temperature variability 
during quiescence is more subtle, with a mean kT of 1 .60 keV 
and a standard deviation of 0.13 keV. 

3.3. Abundances 

In Fig. 5 we plot the relative abundances of the major elements 
with respect to Fe, in solar-abundance units (Asplund et al. 
2009), and as a function of their first ionization potential (FIP). 
In the plot, the data points are progressively shifted along the 
x-axis according to the order in time in which they were mea- 
sured, thus creating a time-like sequence for each element. The 
abundances are expressed relative to the Fe abundance instead 
of the more common absolute values relative to H. The absolute 
Fe abundance with respect to H is 0.32, 0.32, 0.38 and 0.31 so- 
lar in observations 201, 301, 401 and 501Q respectively, with 
a typical error of ±0.03. In the APEC model and the global 
spectral fit approach in the XSPEC software that we use, the 
absolute abundance is constrained by the flux in the Fe lines 
and the overall fit to the continuum. While lines are emitted 
from a limited range of temperautes, continuum flux at all RGS 
wavelengths is contributed from all temperature components. 
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Absolute abundance is therefore sensitive to the underlying tem- 
perature distribution, which is not so well constrained at the high 
temperatures. Conversely, estimates of abundances relative to Fe 
rely on well-measured line ratios that tend to emerge from a lim- 
ited and overlapping range of temperatures. The large number 
of observed Fe ions and lines, whose peak emissivities span a 
wide temperature range from kT — 0.3-1. 5 keV makes these 
measurements significantly less sensitive to the temperature dis- 
tribution in the plasma. The selection of the reference element 
for describing abundances does not affect the general pattern 
of abundance variations. We normalize the abundance ratios to 
those of the solar photospheric composition of Asplund et al. 
(2009) in order to plot them on the same scale and to discuss the 
FIP effect. 

As seen in Fig. 5, El Eri shows the typical inverse FIP ef- 
fect observed in the most active coronae (Brinkman et al. 2001; 
Audard et al. 2003). In particular, the Ne abundance is enriched 
by a factor -3 compared to Fe, and a lesser enrichment is ob- 
served for C, N, and O (all high FIP) over Fe (low FIP), The 
abundance ratio of the low-FIP elements Mg and Si to Fe is ap- 
proximately solar, although with respect to H, they are all sub- 
solar (~0.3). 

Little variation of the abundances with the level of activity is 
observed. The O (high FIP) abundance does seem to be low dur- 
ing the high-states of Obs. 401 and 501F (flare), and high during 
the low-activity state of Obs. 50 IQ. This is reminiscent of the 
chromospheric evaporation effect of solar-composition material 
during flares in an otherwise (quiescent) inverse-FIP coronae, 
which was observed by Nordon & Behar (2007, 2008) in a sam- 
ple of flares. The highly enriched Ne abundance is slightly re- 
duced in the high state of Obs. 401. However, despite the good 
signal, no significant variation in the Ne abundance can be de- 
tected during any other period, not even during the flare, which 
would be expected if chromospheric evaporation takes place. 
The Mg/Fe abundance ratio comes down by -20% during the 
flare, but this can not be a FIP related effect, as both elements 
are low-FIP. 

To better quantify the abundance variations, a single number 
that quantifies the FIP bias can be defined (following Nordon & 
Behar 2008) 

FB = (log(A z /A Zo ) FIP< ioev) - (log(A z /A Z0 ) FIP> ioev) (2) 

where (log (A z /A Z q )) is the error-weighted mean abundance (in 
log) of the low/high FIP elements, relative to a set of reference 
abundances, in this case the solar abundances. We use the tradi- 
tional solar distinction between low FIP (<10 eV) and high FIP 
(>10 eV). A positive FB value indicates a solar-like FIP bias, 
while a negative FB value indicates an IFIP bias. 

The FIP biases in the five observations (20 1, 301, 401, 501 Q, 
50 IF) are shown in Fig. 6. The FB variations over the course of 
the present campaign are not dramatic, but still some trends can 
be identified. The FB in Fig. 6 generally varies with activity, 
namely the FB increases, i.e., tends to (solar) photospheric val- 
ues, as the activity increases, rising from Obs. 201 to Obs. 301, 
and further in Obs. 401, and then decreasing in Obs. 501Q when 
the count rate returns to its lower level. For comparison, see 
the light curves in Fig. 1. This trend breaks down during the 
flare when the activity strongly increases, but the FB decreases. 
At face value, this means that as the (quiescent) coronal ac- 
tivity increases, perhaps more small active regions appear, and 
more high-FIP elements are preferentially brought up into the 
corona. This process is reminiscent of the evaporation of chro- 
mospheric material observed in the largest flares on the most 
active (strongest IFIP) coronae (Nordon & Behar 2007, 2008). 



Fig. 6. FIP bias (FB from Eq. (2)) for each observation and as a func- 
tion of time. Horizontal error bars indicate the exposure time of each 
observation. 

On the other hand, the present flare on El Eri does not follow 
this trend, but the weak abundance variation observed is consis- 
tent with the abundance variations observed in flares on coronae 
with similar quiescent FB values of about -0.2 (see Fig. 5 in 
Nordon & Behar 2008). 

4. Rotation effects 

The four exposures that all cover the same rotation phase al- 
low us to study the dependence of temperatures and abundances 
on phase, by further breaking up each observation into smaller 
segments and co-adding the similar phase segments from each 
observation. If active coronal regions are compact, and sparse 
in longitude, then the rotation would create variations in the ob- 
served X-ray emission properties, as active regions rotate in and 
out of view, over the limb of the star. If the active regions last 
longer than the monitoring time of 6 days (as solar active regions 
do), then stacking data from the same phase, but from different 
orbits, can increase the signal of this rotational modulation. 

During each full observation of -36 ks (Table 1), El Eri ro- 
tates by 72 degrees, i.e., 20% of the orbit. The inclination an- 
gle of the rotation axis to the line of sight of 56° implies only 
about 10% of its surface gradually appears during the observa- 
tion and a different 10% gradually moves out of sight. We divide 
each observation into four quarters of -9 ks each, during each 
of which the El Eri system rotates by 18 degrees. The data are 
then reduced in quarters and the spectra of each instrument are 
combined per phase bin. The two late phase bins of Obs. 501 
are excluded from the analysis, due to the flare (see Fig. 1). 
Subsequently, we perform the same EM and abundance analy- 
sis as before, but for each phase-stacked spectra. 

Figure 7 shows the IEM for each phase bin, Q1-Q4. The 
first three phase bins show the same temperature structure. In 
the last phase bin (Q4), there is an excess of emission at temper- 
atures below -1.5 keV, which slightly increases the total IEM 
of Q4, although it is just consistent with the previous quarters 
to within the errors. This excess emission measure is a result of 
the rise in flux observed in the light curves (Fig. 1) during the 
last quarter of all of the observations, but most conspicuously 
in Obs. 401. A similar rise in flux can be seen also in the last 
quarter of Obs. 501, which is somewhat masked by the flare, 
and therefore excluded from the analysis. Figure 7 can be ex- 
plained by an active region that is not visible at all during the 
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Fig. 7. Cumulatively integrated emission measure distribution (IEM) as 
a function of the four phase bins Q1-Q4. 



Fig. 8. Mean temperature for the four rotation phase bins. Very little 
evolution is observed compared to the days-scale temperature variations 
(Fig. 4). 

two earlier quarters, that comes into view from beyond the limb 
during Q3 or Q4 (i.e., after ~40 degrees of rotation), and that 
gradually brightens during the four days after the beginning of 
the XMM-Newton monitoring, i.e., it becomes particularly bright 
during Obs. 401 and Obs. 501. 

In Fig. 8 we plot the mean temperature as a function of phase. 
The value of kT e g- is very stable between the observed rotation 
phases, despite the added emission measure. In fact, it shows 
much less scatter than the scatter from one observation to the 
other over several days (Fig. 4). The lack of temeprature vari- 
ation is somewhat expected, since each phase represent a time- 
average of 4-6 days, and variations on shorter time scales tend to 
average out. Note that here too, the two last quarters of Obs. 501 
are excluded from the analysis due to the flare. The emerging ac- 
tive region (presumably a magnetic loop) reaches the mean tem- 
perature of (kT ~ 1.5 keV), which is very similar to the mean 
temperature of the coronal regions on the rest of the El Eri ob- 
served area. 

The mean abundances for each rotation phase bin are plotted 
in Fig. 9. Again, the data points are progressively shifted along 
the horizontal axis from left to right, in order to give an idea of 
the evolution with phase. The variation in abundance with phase 
for each element is very small, and statistically insignificant 



FIP [eV] 


Fig. 9. Coronal abundances relative to Fe and expressed relative to the 
solar X/Fe abundance ratio (Asplund et al. 2009). Points from different 
phase bins Q1-Q4 are progressively shifted to the right along the hor- 
izontal axis to enhance visibility and to emphasize the variability with 
phase. 



Fig. 10. FIP bias ( FB ) as a function of rotation phase averaged over all 
available observations in four iso-phase intervals. Horizontal error bars 
reflect the ~9 ks duration in each bin. 

between consecutive bins. However, Mg, O and Ne, the three 
elements with the best signal do appear to vary monotonously 
with phase. The abundance ratio of Mg/Fe decreases while that 
of O/Fe and Ne/Fe increases. The other elements tend to show 
a similar monotonic variation, although with much larger uncer- 
tainties. While the trend is of low signal in each element, the re- 
peated pattern in several elements makes it more believable, and 
less likely due to noise. Moreover, since the trend of the low-FIP 
Mg is opposite to that of the high-FIP Ne and O, we suspect it 
is due to a FIP related effect and not due to only Fe changing 
with phase. 

As in Sect. 3.3, we wish to enhance the individual abundance 
signals as to reflect potential FIP related effects by computing 
the FIP bias from Eq. (2) for each phase-stacked spectra. The 
FIP bias as a function of phase is plotted in Fig. 10. The FB 
shows a mild, yet monotonous decrease with phase. This can 
be interpreted as the emerging active region contributing an in- 
verse FIP bias. In other words, the active region (Q3, Q4) has 
slightly more high-FIP elements than the background corona 
(Ql, Q2). Note that this time-averaged rotation effect is different 
from the time-evolution observed over days (Fig. 6), in which the 
FB increases with increased activity. Here, the increased activity 
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coming from a specific region on the star is the likely cause of 
the decrease in FB. 

The weakness of the observed signal is not unexpected, as 
the fraction of the El Eri surface rotating in and out of our view 
is not very large, ~5% between quarters. Say we measure during 
rotation phase Q4, in which the active region is visible, a mean 
abundance A (4) which deviates from the mean abundance dur- 
ing phase Q1 A z } (in which the active region is not visible) by 
a ~ 0.1, i.e. A (4) = (1 + a)A z \ The total count rate is increasing 
with phase, indicating that this effect is likely due to a different 
abundance in the active region which is rotating into view, for 
which A™ 810 " = (1 + fl)A z \ Since the contributions to the spec- 
trum of both the background and the active-region abundances 
are weighted by their respective emission measures, and we at- 
tribute all the excess emission measure during Q4 over Q1 to the 
active region, then one can easily show that: 


where EM* 4 * and EM* 1 * are the total emission measures in 
Q4 and Q1 respectively. According to Fig. 7, the EM ratio is 
EM*' 4 ' 1 / EM 1 ~ 1.15. Therefore, the actual abundance variation 
in the active region by comparison with the spatial-mean back- 
ground corona is at the ~ 0.75 level, or a 75% stronger IFIP ef- 
fect than the mean corona. Admittedly, this conclusion for El Eri 
is based on a marginal signal and is, thus, associated with ap- 
preciable uncertainty, but it does point to significant abundance 
variations in a compact active region with respect to the mean 
corona. 

5. Discussion and conclusions 

The light curves and abundance patterns in the corona of El Eri 
show an active region appearing over the limb as the corona 
rotates with the star. Although on time scales of days, the in- 
crease of the activity level seems to lead to an increase of the FB 
(i.e., less IFIP bias, abundances more akin to the photosphere), 
the appearance of the active region with phase features an in- 
crease in high-FIP elements and a lowered FB (i.e. a tendency 
to a stronger IFIP bias). Two different effects are at play: one is 
the evolution of the emission measure and coronal abundances 
on a time scale of days between exposures, spatially integrated 
and averaged over the entire visible surface of the star. The other 
is an already-existing, compact and bright active-region that per- 
sists over the six days duration of the monitoring campaign and 
rotates into view towards the end of each exposure. Thus, the 
active regions is manifested as variations with phase instead of 
time. 

One cannot tell whether the days-scale variations in count 
rates, emission measure, and abundances are due to the evo- 
lution of existing coronal loops, or due to the appearance of 
new coronal loops and disappearance of old loops. However, 
the correlation in which an increase in emission measure and 
count rates corresponds to a decrease in the IFIP bias (increase 
in FB, closer to photospheric composition) offers an interesting 
interpretation. If the brightening is due to the emergence of new 
loops, and those appear initially with nearly photospheric com- 
position, then the mean coronal abundances will tend towards 
photospheric composition during these times. Conversely, dur- 
ing times in which no new loops emerge and existing loops 
age, evolve in composition and perhaps some fade, the mean 
abundances will tend to resemble the composition of the older, 
evolved loops, i.e. different than photospheric. 


Indeed, on the sun, it has been reported that solar coronal 
loops evolve on time scales of days (e.g., Widing & Feldman 
2001). Newly emerging structures tend to have photospheric 
composition, but develop the solar FIP bias within a few days. 
Moreover, on the sun, the typical enhancement of low FIP abun- 
dances relative to photospheric values by a factor of 4 is reached 
within 2-3 days and continues to evolve further, reaching a fac- 
tor of 8 or more. The analogous interpretation and one in agree- 
ment with the present observations is that new coronal loops on 
El Eri emerge with photospheric composition and develop an 
IFIP bias over time scales of days. The spatial mean IFIP bias is 
determined by mean age of the loops. Thus, periods of gradual 
brightening of the corona are characterized by the emergence of 
young loops that reduce the global IFIP bias, while periods of 
gradual dimming are characterized by aging of existing loops 
that increase the global IFIP bias. Note that the negative correla- 
tion between emission measure and IFIP bias as described above 
relates to the days scale evolution of a given stellar corona and is 
not at odds with the general positive correlation found between 
activity level and tendency to IFIP bias in stellar populations. 

Rotational modulation in X-rays and EUV has been re- 
ported for HR 1099 (Agrawal & Vaidya 1988; Drake et al. 1994; 
Audard et al. 2001), from light curves, but phase dependent spec- 
troscopy was not possible, due to the constant flaring of that 
source, which makes it difficult to disentangle between rotational 
modulation of active regions and their intrinsic activity variabil- 
ity. A different, yet notable, phase dependent measurement is 
that of line profiles by Brickhouse et al. (2001), which enabled 
latitude determination of active coronal regions on the stellar bi- 
nary 44i Boo. The present measurement, however, is the first to 
measure phase dependent abundances. 

The active region that appears towards the end of each ex- 
posure persists throughout the six days of the observing cam- 
paign. It therefore likely contains loops that are at least a few 
days old. In the scenario which we proposed above, these loops 
will have already evolved (over days, possibly starting even be- 
fore the beginning of the campaign) and their abundances devel- 
oped a strong IFIP bias, stronger than the coronal mean bias that 
includes younger loops. Thus, as this region comes into view, the 
spatially averaged abundances tend towards a stronger IFIP bias 
(Fig. 10). 

After showing that large flares tend to counter the coronal 
FIP effect (Nordon & Behar 2007, 2008) - the present work 
promotes the notion that the abundances in individual coro- 
nal active regions, differ from the mean coronal abundances as 
well. It is not clear whether large flares show coronal abun- 
dances that are closer to photospheric composition than the 
quiescent corona because of added material from the photo- 
sphere/chromosphere (chromospheric evaporation) or because 
these flares happen preferentially in young active regions where 
the abundances are closer to photospheric than the global mean. 
The abundances measured during the flare observed in this cam- 
paign show an increase IFIP bias (more negative FB, Fig. 6) in- 
stead of the typical tendency to bring abundances towards pho- 
tospheric composition. If the flare is related to the limb active 
region that is strongly IFIP biased, it is possible that the mean 
abundances, weighted between the flaring active region, evapo- 
rated chromospheric material and surrounding quiescent corona 
end up being dominated by the abundances of the active region. 

While the proposed scenario in which young loops emerge 
with photospheric-like composition and then develop over time- 
scale of days to an IFIP bias can explain our observations, it does 
not explain why the abundances evolve an IFIP effect instead of 
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Fig. A.l. RGS spectrum (RGS1 and RGS2 combined) of Obs. ID 060317201. Plotted model is a simultaneous fit to all instruments (see EPIC 
spectra in Fig. A.2). 



Fig. A.2. EPIC-pn and EPIC-MOS1 spectra of Obs. ID 060317201. Plotted model is a simultaneous fit to all instruments (see RGS spectrum in 
Fig. A.l). 


a solar-like FIP effect. For a theory that attempts to resolve this 
issue see Laming (2004, 2009, 2012). 
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Appendix A: Spectra 

Sample spectra and best-fitted model are presented in this sec- 
tion. The X-ray spectra obtained from our first observation 
ID 060317201. as measured by all XMM-Newton instruments 
used in the analysis are shown in Fig. A.l. This observation 
has the lowest count rates (see Fig. 1). Consequently, the spec- 
tra of the three other observations have slightly better signal-to- 
noise ratios, but are overall similar to those shown in Fig. A. 1 . 
In Sect. 4, each observation is divided into four segments, and 
the four iso-phase segments are then stacked, resulting in spec- 
tra that are, again, with numbers of photons and signal-to-noise 
ratios that are comparable to those of the individual observations. 

Appendix B: IEM and ( kT ) error estimates 

The EMD solutions to coronal emission tend to be highly degen- 
erate as emission measure can be traded between different, but 
close temperatures with little effect on the resulting spectrum. 
Thus, when representing the EMD with isothermal components, 


the uncertainties in the normalizations of the various compo- 
nents are highly correlated, especially when many components 
are used. The correlation tends to be negative and increases as 
the components are closer in temperature, i.e., adjacent temper- 
atures can compensate for each other. This results in large errors 
on each individual temperature component, which can be seen 
in Fig. 2. On the other hand, the exchange of emission measure 
between close-temperature components has only a small effect 
on the cumulative, or integrated emission measure (IEM). Still, 
the spectral line emissivities are temperature dependent and thus 
the IEM is not perfectly conserved. In addition, there are the un- 
correlated uncertainties associated with the limited statistics of 
the data. To estimate the highly nonlinear contributions of these 
effects on the uncertainties in the IEM, we allude a Monte-Carlo 
simulation. 

The fitting algorithm in XSPEC gives, along with the best 
fit model, the covariance matrix of second derivatives calculated 
at the best-fit solution. Here, we use this covariance matrix to 
redraw a random, new set of model EMD values on the given 
temperature grid (see Fig. 3) that is close in its statistical prop- 
erties to the original best-fit model, in the sense that it produces 
spectral fits to the data that are almost as good as the best-fit 
model. To generate a random set of numbers corresponding to 
a correlation matrix R, one can decompose R into R = U T U. 
If X is a vector of uncorrelated random values, X c = XU is a 
set of correlated random values corresponding to R. The process 
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of generating random EMDs is repeated numerous times, and 
for each EMD solution the corresponding IEM(kT) is calculated 
at each temperature point on the grid. Thus, we generate a large 
sample of IEM(kT) values at each temperature point, from which 
we adopt the central 68% as the equivalent of the 1-cr uncertainty 
interval that is, then, plotted in Figs. 3 and 7. For each generated 
random EMD we also calculate the emission-measure weighted 
mean temperature ( kT ). The central 68% of the (kT) distribu- 
tion is adopted as the 1-cr uncertainty interval on (kT), plotted 
in Figs. 4 and 8. 
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